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Abstract — We study spanning diverging forests of a digraph and related matrices. It is shown 
that the normalized matrix of out forests of a digraph coincides with the transition matrix in a 
specific observation model for Markov chains related to the digraph. Expressions are given for 
the Moore-Penrose generalized inverse and the group inverse of the Kirchhoff matrix. These 
expressions involve the matrix of maximum out forests of the digraph. Every matrix of out 
forests with a fixed number of arcs and the normalized matrix of out forests are represented as 
polynomials of the Kirchhoff matrix; with the help of these identities, new proofs are given for 
the matrix-forest theorem and some other statements. A connection is specified between the 
forest dimension of a digraph and the degree of an annihilating polynomial for the Kirchhoff 
matrix. Some accessibility measures for digraph vertices are considered. These are based on 
the enumeration of spanning forests. 

1. INTRODUCTION 

Directed graphs provide a simple and universal tool to model connection structures. It is not 
accidental that the first systematic monograph in the theory of digraphs [1] was titled "Structural 
Models: An Introduction to the Theory of Directed Graphs." Digraphs frequently serve to model 
processes that can proceed in the direction of arcs. Physical transference, service, control, trans- 
mission of influences, ideas, innovations, and diseases are examples of such processes. If a process 
can start from a number of vertices and ends with the inclusion of all vertices, then the process 
can be modelled by the family of out forests (i.e., spanning diverging forests) of the digraph. The 
enumeration of all out forests allows one to determine the typical roles of the vertices in the process: 
one vertex is a typical starting point, another vertex is a typical intermediate point, some vertex 
is a typical terminating point of the process, etc. If an initial (weighted) digraph imposes some 
measure on the said processes, then the "role profile" of each vertex can be expressed numerically. 
Moreover, an exact answer can be given to the following important question: how likely is it that 
the process initiated at vertex j arrives at vertex i. It is not surprising that out forests of a digraph 
turn out to be closely related with Markov chains realizable on the digraph. 

The study of out forests has been started in [2] . Generally, they were given less attention in the 
literature, than that given to spanning diverging trees (out arborescences), which exist only for a 
narrow class of digraphs. We mention in this connection [3-11], where still undirected forests were 
considered in most cases. The maximum out forests (i.e., out forests with the greatest possible 
number of arcs) of a digraph were studied in [12,13]. It was established that the normalized matrix 
of such forests coincides with the matrix of limiting probabilities of every Markov chain related to 
the given digraph. Some results on spanning forests of directed and undirected multigraphs were 
given in [14, 15]. 

1 This work was supported by the Russian Foundation for Basic Research and INTAS. 
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In this paper, we study the normalized matrix of out forests (which has been also termed the 
matrix of relative forest accessibilities and the matrix of forest proximities) and the matrices of 
forests with fixed numbers of arcs. 

2. NOTATION AND SOME EARLIER RESULTS 

In the terminology, we mainly follow [1,16]. Suppose that J 1 is a weighted digraph without 
loops, V(r) = {1, . . . , n} (n > 1) is its set of vertices, and E{T) its set of arcs. The weights of all 
arcs are supposed to be strictly positive. A subgraph of a digraph U is a digraph whose vertices 
and arcs belong to the sets of vertices and arcs of r; the weights of subgraph's arcs are the same as 
in r. A restriction of r to V' C V(r) is a digraph whose arc set contains all the arcs in E(r) that 
have both incident vertices in V' . A spanning subgraph of r is a subgraph with vertex set U(U). 
The indegree id(iw) of vertex w is the number of arcs that come to w, outdegree od(u>) of vertex w 
is the number of arcs that come from w. A vertex w will be called undominated if id(ty)=0 and 
dominated if id(w) > 1. A vertex w is isolated if r contains no arcs incident to w. 

A route in a digraph is an alternating sequence of vertices and arcs wo,e±, w±, . . . , e^, Wk with 
every arc ei being (wi-i,Wi). If every arc is either (wi-i,Wi) or (wi, Wi-±), then the sequence is 
called a semiroute. A path in a digraph is a route all whose vertices are different. A circuit is a 
route with w = Wk, the other vertices being distinct and different from u>o. A vertex w is reachable 
from a vertex z in r if w = z or F contains a path from z to w. A semicircuit is an alternating 
sequence of distinct vertices and arcs, wo,ei,wi, . . . ,ek,wo, where every arc ei is either (wi-\,Wi) 
or (wi,Wi-i) and all vertices wq, . . . , Wk-i are different. The restriction of r to any maximal subset 
of vertices connected by semiroutes is called a weak component of r. Let E = (£ij) be the matrix 
of arc weights. Its entry is zero if and only if there is no arc from vertex i to vertex j in r. If 
r' is a subgraph of r, then the weight of r', e(r'), is the product of the weights of all its arcs; if 
r' does not contain arcs, then e(r') = 1. The weight of a nonempty set of digraphs Q is defined as 
follows: 

e{Q) = ]T e(H); 
Heg 

the weight of the empty set is 0. 

The Kirchhoff matrix [17] of a weighted digraph r is the n x n-matrix L = L(r) = (^ -) with 

elements £^ = —En when j ^ i and i u = — Yl ^iki h J ' = 1> • • • > n - 

k^i 

A diverging tree is a digraph without semicircuits that has a vertex (called the root) from which 
every vertex is reachable. The indegree of every non-root vertex of a diverging tree is 1. If w is the 
root, then id(iu) = 0. A converging tree is a digraph without semicircuits that has a vertex (called 
the sink) reachable from every vertex. 

A diverging forest {converging forest) is a digraph without circuits such that \d(w) < 1 (respec- 
tively, od(w) < 1) for every vertex w. An out forest (in forest) of a digraph r is any its spanning 
diverging (respectively, converging) forest. 

The weak components of diverging forests (converging forests) are diverging trees (respectively, 
converging trees). 

Definition 1. An out forest F of a digraph r is called a maximum out forest of r if r has no 
out forest with a greater number of arcs than in F. An in forest F of a digraph r is a maximum 
in forest of r if r has no in forest with a greater number of arcs than in F. 

Obviously, every maximum out forest of r has the minimum possible number of weak compo- 
nents (out trees) ; this number will be called the out forest dimension of the digraph and denoted 
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by v. The number of arcs in any maximum out forest is obviously n — v. The number of weak 
components of every maximum in forest will be called the in forest dimension of the digraph and 
denoted by v'. Obviously, for every digraph, v,v' E {1, . . . , n}. 

If a digraph i~i is obtained from r by the reversal of all arcs, then the out forests in r naturally 
correspond to the in forests in i~i and vice versa. Therefore, the out forest dimension and in forest 
dimension of r are respectively equal to the in forest dimension and out forest dimension of I\. 

The following proposition states that the dimensions v and v' of a digraph are not connected, 
except for the case where v = n and v' = n. 

Proposition 1. 1. Let k,k' G {l,...,n — 1}. Then there exists a digraph on n vertices such 
that v = k and v' = k' . 

2. For every digraph r on n vertices, v = n<$v' = n<3- E{r) = 0. 

The proofs are given in the Appendix. 

Throughout the paper, we mainly deal with diverging forests. However, all the results have 
counterparts formulated in terms of converging forests. Simple properties of out forests have been 
studied in [13] (Section 3). We do not cite them here and only confine ourselves to the following 

Proposition 2. If i and j belong to different trees in a maximum out forest F of a digraph r, 
and j is a root in F, then r contains no paths from i to j. 

Let us adduce some definitions and results from [13] which are frequently used below. 

Definition 2. A nonempty subset of vertices K C V(r) of digraph r is an undominated knot 2 
in r iff all the vertices that belong to K are mutually reachable and there are no arcs (vjj,Wi) such 
that Wj G V(r) \ K and Wi G K. 

u 

Suppose that K = [_] Ki, where K±, . . . , K u are all the undominated knots of r, and Kf is the 

i=i 

set of all vertices reachable from Ki and unreachable from the other undominated knots. For any 
undominated knot K of r, denote by Fx the restriction of F to K and by F-k the subgraph with 
vertex set V(r) and arc set E(r) \ E{Fk)- For a fixed K, T will designate the set of all spanning 
diverging trees of Fk and V will be the set of all maximum out forests of F-k- By T k , k G K, we 
denote the subset of T consisting of all trees that diverge from k, and by j? K ~* t j % £ V(r), the set 
of all maximum out forests of F_k such that i is reachable from some vertex that belongs to K in 
these forests. 

By !F{r) = T and Tk{E) = Tk we denote the set of all out forests of F and the set of all out 
forests of F with k arcs, respectively; will designate the set of all out forests with k arcs where 

j belongs to a tree diverging from i. 

Definition 3. The matrix J = (Jy) = a~ 1 Q n - v , where a = e(T n - v ), Q n -v = ( a ij) = 
(£(.Fn^>))> W ^ k e cane d the normalized matrix of maximum out forests of a digraph. 

Theorem 1. Suppose that r is an arbitrary digraph and K is an undominated knot in r . Then 
the following statements are true. 

_ n 

1. J is a stochastic matrix: Jy > 0, J2 Jik = lj i,j = l,---,n- 

fc=i 

2. Jij 7^ O (j G K and i is reachable from j in T). 
2 In [2], undominated knots are called W-bases. 
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3. Suppose that j G K. For any i € V{T), Jij = e(T j )e(V K ^ i )/e(F n - v ). Furthermore, if 
ie K+, then = Jjj = e(T j )/e{T). 

4- J2 Jjj = 1- I n particular, if j is an undominated vertex, then Jjj = 1. 
jeK 

5- Ifji,j2 £ K, then J.j 2 = (e{T^ 2 ) / e{T^ 1 )) J.j 1 , i.e., the j\ and j2 columns of J are proportional. 
Theorem 2. For every weighted digraph, J is idempotent: J 2 = J . 

Theorem 3. For every weighted digraph, LJ = JL = 0. 

Theorem 4 (a parametric version of the matrix- forest theorem). For any weighted multidigraph 
r with positive weights of arcs and any r > 0, there exists the matrix Q(t) = (I + rL(r))^ 1 and 

-, n—v 

where 

s(r) = J2r k e(T k ), Q k = (qf J ), = e{J^), k = 0,...,n-v, i,j = l,...,n. (2) 

k=0 



Definition 4. The matrix Q k , k = 0, . . . , n — v , will be called the matrix of out forests of r 
with k arcs. 

Theorem 4 represents (/ + rL) _1 via the matrices of out forests with various numbers of arcs. 

Definition 5. The matrices Q(r) = (I + tL)^ 1 , t > 0, will be called the normalized matrices 
of out forests of a digraph. 

In [14], the matrices Q(t) = (I + tL) _1 were referred to as the matrices of relative forest 
accessibilities of a digraph. In Section 4, Q(r) are expressed as polynomials of L (Corollary from 
Theorem 8). 

Theorem 5. For every weighted digraph r, lim Q(t) = lim (I + rLH = J . 

T — *O0 T—*CO 

3. MATRICES OF OUT FORESTS AND TRANSITION PROBABILITIES OF MARKOV 

CHAINS 

It has been shown in [13] that the matrix of Cesaro limiting probabilities of a Markov chain 
coincides with the normalized matrix J of maximum out forests of any digraph related to this 
Markov chain. Now we give a Markov chain interpretation for the normalized matrices of out 
forests Q(t) with any r > 0. 

Definition 6 [13]. A homogeneous Markov chain with set of states {1, . . . ,n} and transition 
probability matrix P is related to a weighted digraph r iff there exists a/0 such that 

P = I-aL(r). (3) 
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Let r be a weighted digraph. Consider an arbitrary Markov chain related to r and the following 
observation model. 

The geometric model of random observation. Suppose that a Bernoulli trial is performed 
at the point of time t = with success probability q (0 < q < 1). In case of success, t = becomes 
the epoch of observation. Otherwise, Bernoulli trials are performed at t = 1, 2, . . . — to the point of 
the first success. This point becomes the epoch of observation. 

This model determines a discrete probability distribution p{k) of the epoch of observation on the 
set {0, 1, 2, . . .}. This is obviously the geometric distribution (which gives the name of the model) 
with parameter q: 

p(k)=q(l-q) k , fc=l,2,... (4) 

Consider Markov chain multistep transitions in a random number of steps: from the initial state 
at t = to the state at the random epoch of observation distributed geometrically with parameter q. 

Suppose that P(a, q) = (p^ (a, q)j is the matrix of unconditional probabilities for such multistep 
transitions: from the initial state to the state at the epoch of observation. 

Theorem 6. For any weighted digraph, any r > and any Markov chains related to the weighted 
digraph, 

Q{r)=P{a,q) 

holds, where 

q=(T/a + l)-\ (5) 

Theorem 6 provides an interpretation for the normalized matrix Q(t) of out forests in terms 
of Markov chain transition probabilities. Conversely, for any Markov chain, the transition proba- 
bilities in the geometric observation model can be interpreted in terms of diverging forests of the 
corresponding digraphs. 

The following corollary stresses the arbitrariness of Markov chains in Theorem 6. 

Corollary 1 from Theorem 6. For every Markov chain, every success probability q e]0, 1[ 
in the geometric observation model, and every digraph related to the Markov chain, 

P(a,q)=Q(r) 

holds, where r = — l)a. 

Corollary 2 from Theorem 6. 

-. fc-i 

lim P(a, q) = J= lim - V P k . (6) 



By Corollary 2 from Theorem 6, at a vanishingly small success probability q, the transition 
probabilities in the geometric observation model are given by the matrix J of maximum out forests 
of any weighted digraph to which this chain is related. 
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4. REPRESENTATIONS OF FOREST MATRICES VIA THE KIRCHHOFF MATRIX AND 

THEIR CONSEQUENCES 

In this section, we represent the matrices Q k of out forests with k arcs as polynomials of the 
Kirchhoff matrix L (Theorem 7). This allows one to obtain alternative proofs of Theorems 2-4 and 
to represent the matrix Q(t) = (I + tL)^ 1 as a polynomial of L (Theorem 8). Proposition 4 gives 
an easy way to calculate Q k , k = 1, . . . , n — v , and J. 

By a k we denote the total weight of all out forests of r with k arcs: a k = e(F k ), k = 0, . . . , n — v. 

Proposition 3. For any weighted digraph and any k = 0, . . . , n — v, 

Qk+i = o- k+1 I - LQ k . (7) 

Observe that since the weight of the empty set is 0, we have Q n - v+ \ = and a n _ v+1 = 0. 
Taking the traces on the left-hand side and the right-hand side of (7) and using the fact that 

MQfc) = (n - k)a k , k = 0, . . . , n - v + 1 (8) 

(because every out forest with k arcs has n — k roots), we deduce 

tr(LQ k ) 

°k+i = k - 1 , k = 0,...,n-v. (9) 

Substituting (9) in (7) provides 
Proposition 4. For every weighted digraph, 

Q k+ i = t -^f^I-LQ k , k = 0,...,n-v. (10) 

Identity (10) enables one to recursively determine the matrices Q k , k = 0, . . . ,n — v, and J, 
starting with Qq = I. Note that this procedure essentially coincides with Faddeev's algorithm [19] 
for the computation of the characteristic polynomial as applied to L. Thus, the matrices involved 
in Faddeev's method are precisely Q k - 

From Proposition 3, it follows 

Theorem 7. For any weighted digraph and any k = 0, . . . , n — v, 

k 

Q k = J2a k _ i (-L) i . (11) 

i=0 

Corollary 1 from Theorem 7. For every weighted digraph, matrices Q k , k = 0, . . . , n — v , 

commute with all matrices with which L commutes, in particular, with L, J, Q{t), and each other. 



Lemma 1. For any k = 0, . . . , n — v, every row sum of LQ k is 0. 
From Proposition 3 and Lemma 1, it follows 

Proposition 5. The matrices LQ k , k = 0, . . . , n—v, are the Kirchhoff' matrices of some weighted 
digraphs. 
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Corollary 2 from Theorem 7. For any weighted digraph, LQ n ^ v = Q n - V L = 0. 

In view of Definition 3, this corollary is equivalent to Theorem 3. Thus, we get a new proof of 
this theorem. 

Consider the matrices 

■h = Ok 1 Qk, k = 0, . . . ,n - v. (12) 
In particular, Jq = I and J n - V = J- 

Making use of the last corollary, we obtain 

Corollary 3 from Theorem 7. For any k £ {1, . . . ,n — v}, J^J = JJk = J- In particular, 
J n -v J = J = J ■ Moreover, Q{t)J = JQ{t) = J for every r > 0. 

This corollary provides a new proof of Theorem 2. 

By virtue of Proposition 3 and Corollary 1 from Theorem 7, the matrices Jk are connected as 
follows: 

Jk+i = I -^JkL, fc = 0,...,n-u-l, (13) 

a k+i 

and, by Lemma 1, each their row is unity. The entries of Jk are nonnegative by definition, thus, 
we obtain 

Proposition 6. For every weighted digraph r, matrices Jk, k = 0, . . . , n — v, are stochastic. 

Completing Proposition 3 with the obvious equality Qo = I = a 1 gives 

Qo = o- i, 

Qi + LQ = cr 1 1, ^ 

Qn—v ~\~ LQn—v—1 = 0~n—v ^ • 

Add up these equations and, using Corollary 2 from Theorem 7, substitute (J + L)Q n _ v for 

Qn—v'- 

n—v 

(I + L)Q + (I + L)Qi + ... + (/ + L)Q n _ v = (J2 <r k )l. 

k=0 

n—v 

Making use of the nonsingularity of I + L (Theorem 4) and the notation s = e{F) = J2 a h' we 

k=0 

obtain 

Y / Q k = s(I + L)-\ (15) 

k=0 

which provides 

Corollary from Proposition 3. For any weighted digraph, 

Q(l) = (/ + L)- 1 = s - in j2Qk. 

k=0 

This statement coincides with the matrix-forest theorem for digraphs [14] and with Theorem 4 
in the case of r = 1. Accordingly, we obtain a new proof of the matrix-forest theorem. 

n—v 

By means of Theorem 7, the matrices Q(t) = (I + tL)~ = s~ l T k Qk including Q(l) = 

k=0 

(I + L) 1 can be represented as polynomials of L. 
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Theorem 8. For any weighted digraph, 

n—v 

Q(l) = (/ + L)- 1 = s- 1 ]T s m (-L)\ (16) 

i=0 

k 

where s k = J2 °~j ^ the total weight of out forests of r with at most k arcs, k = 0, . . . , n — v. 

3=0 

Corollary from Theorem 8. For any weighted digraph and any r > 0, 

n—v 

Q{r) = (I + tL)- 1 = s-\r) ]T s n _ v _ t (r)(-rLy , (17) 

i=0 

k 

where s k (r) = J2 t 3 °~ j-, k = 0, . . . ,n — v. 

3=0 

Note that s(I+L)~ l is the adjugate (the transposed matrix of cofactors) of I+L; s(t) (I+tL)~ 1 
is the same for I+tL. Theorem 8 and the above corollary provide representations for these matrices 
as polynomials of L: 

n—v 

»(j+i)" 1 = Ev«-iN i 1 

i=Q 
n—v 

s(r) (I + tL)- 1 = £ s n _ v _ t (r) {-tL)\ (18) 

i=0 

Remark 1. Since LQk the is Kirchhoff matrix of some weighted digraph (Proposition 5), all 
its principal minors are nonnegative (by Theorem 6 in [2]). Therefore, all LQ^ are singular M- 
matrices (see, e.g., item (Al) of Theorem 4.6 in [20]). Alternatively, this can be concluded from the 
nonnegativity of the real parts of the eigenvalues (see Proposition 9 below) and the nonpositivity 
of off-diagonal elements of L (item (F12) of Theorem 4.6 in [20]). It follows from the representation 
o~ k+1 I — Qk+i = LQk (Proposition 3) of the singular M-matrix LQ^ that o~ k+1 = p(Qk+i), i- e -> 
o~ k+l is the spectral radius of Qk+i, k = 0, . . . ,n — v — 1. This also follows from Proposition 6 
(see (12)). 

5. ON SOME LINEAR TRANSFORMATIONS RELATED TO DIGRAPHS 

For a matrix A G R nxn , by A we denote the linear transformation A : W 1 — > R n induced by A 
with respect to the standard basis of M n : A(x) = Ax. 7£(A) and j\f(A) will designate the range 
and the null space of A, respectively. 

As has been seen in [13], the dimensions of 1Z(J) and 1Z(L T ) are v and n — v, respectively. 
Furthermore, 7£(L T ) n TZ(3) = {0} and, since the dimensions of 1Z(L T ) and 1Z(J) sum to n, M. n 
decomposes to the direct sum of TZ(L T ) and 1Z(J): 

R n = 1Z(L T )+TZ(3). (19) 

Since LJ = (Theorem 3), we get j\f(L) = K(J) and M(J T ) = H(L T ), thus, the sum (19) is 
orthogonal. 

Similarly, in view of JL = 0, the orthogonal decomposition 

R n = TZ(L)+TZ{3 T ) 
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holds along with K(L) n K(J T ) = {0}, M{3) = K(L), and N(L T ) = K(J T ). 

In accordance with (19), every vector u G M n is uniquely represented as u = u x +u 2 , where u 1 G 
K(L T ) = M(J T ) and u 2 G ft(J) = Af(L). For every u / 0, we have (L + J T )u = (L + J 7 )^ + 
(L + J T )u 2 = Lu 1 +J T u 2 . If Lu 1 + J T u 2 = 0, then, since TZ(L) f) TZ(J T ) = {0}, we have Lu x = 
Ju 2 = 0, whence, by N(L) nJV(J ) = {0}, u 1 = u 2 = results. Therefore, the dimension of the 
range (rank) of Z = L + J T is n. Thus, we obtain 

Theorem 9. For any weighted digraph r, the matrix Z = L + J T is nonsingular. 

We will also need the nonsingularity of L + J. 

Theorem 10. For any weighted digraph r, the matrix L + J is nonsingular. 

Corollary from Theorem 10. For any weighted digraph and any a / 0, the matrix L + a J 
is nonsingular. 

~ 2 — — 

It follows from J = J (Theorem 2) that every nonzero columns of J is an eigenvector of J 
associated with the eigenvalue 1. Hence, for any u G 7£(J), J(u) = u holds, therefore, TZ(J) is 
exactly the subspace of fixed vectors of J. 

6. THE MOORE-PENROSE AND GROUP INVERSES OF THE KIRCHHOFF MATRIX 

In this section, we obtain some expressions for the Moore-Penrose generalized inverse and the 
group inverse of the Kirchhoff matrix L. The Moore-Penrose generalized inverse of a rectangular 
complex matrix A is the unique matrix X such that 

(1) AX A = A, 

(2) XAX = X, 

(3) (AX)* = AX, 

(4) (XA)* = XA, 

where (AX)* and (XA)* are the conjugate transposes (Hermitian adjoints) of AX and XA, re- 
spectively. 

For any matrix A, the Moore-Penrose generalized inverse, A + , does exist and is unique. If A is 
nonsingular, then A + coincides with A^ 1 . 

The Moore-Penrose inverses are of theoretical and practical interest. The latter is because A + 
provides the normal pseudosolution of the inconsistent equation ^4x = b: it is x = The 
normal pseudosolution is a vector of the minimum length that minimizes the length of Ax. — b (the 
minimum norm least-squares solution). As applied to Laplacian matrices, such solutions, among 
others, were considered for some preference aggregation problems (more specifically, estimation 
from paired comparisons) [23] , in constructing geometrical representations for systems modelled by 
graphs [24], in the analysis of social networks, and cluster analysis. 

The group inverses are no less important (see, e.g., [22]). A matrix X is the group inverse of 
a square matrix A, if X satisfies the conditions (1) and (2) in the definition of Moore-Penrose 
generalized inverse and also 

(5) AX = XA. 

The group inverse of A is denoted by A*. Generally, group inverses need not exist, but if such 
a matrix exists, then it is unique, but A + = A* is not necessary. 

If L is symmetric (in particular, this is the case for symmetric digraphs, which can be identified 
with undirected graphs), then the matrix (L + a J) -1 — a -1 J (with any a > 0) is [15] the Moore- 
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Penrose generalized inverse and the group inverse of L. Moreover, the latter is true for every 
digraph. 

Theorem 11. For every weighted digraph and any 

L* = (L + aJ)~ 1 -a' 1 J (20) 

and 

L*L = LL* = 1 - J. 

As well as in the case of undirected graphs, L* = {if-) can be obtained via a passage to the 
limit. 

Proposition 7. For every weighted digraph, 

L*=Mmr(0(r)-j). 

We now express L* in terms of the normalized matrices of out forests J n -v-i and J n - V = J 
(see (12)). The following proposition is an analogue of Theorem 3 in [15]. 

Proposition 8. For every weighted digraph, 

L* = ( J n - v -\ — J) . 

°~n—v 

Because of the nonsymmetry of / — J = L*L = LL*, L* is not generally the Moore- Penrose 
generalized inverse of L for digraphs. To obtain an explicit formula for L + , consider the matrix 
Z = L + J T which, by Theorem 9, is nonsingular. Using the identity L J = (Theorem 3), we 
obtain 

{z T y l z- 1 = {zz T y l = (j t j+ll t )-\ 

Lemma 2. For every weighted digraph, (ZZ T )^ 1 commutes with LL T and J 1 J '. 

Matrices LL T , J T J, and (ZZ T )~ 1 are symmetric. The product of two symmetric matrices is 
symmetric iff they are commuting [25] . This implies the following corollary. 

Corollary from Lemma 2. For every weighted digraph, the matrices LL T (ZZ T )^ 1 and 
J T J{ZZ T )~ l are symmetric. 

These facts are useful for the proof of the following theorem. 

Theorem 12. For every weighted digraph, the matrix L T (ZZ T )~ 1 = L T (J T J +LL T )~ l is the 
Moore-Penrose generalized inverse of L. 

7. ON THE GERSGORIN REGION AND ANNIHILATING POLYNOMIALS FOR THE 

KIRCHHOFF MATRIX 

By the Gersgorin theorem (see, e.g., [25]), the eigenvalues of a matrix A belong to the union 
G(A) of n discs: 

G(A) = {j{zeC\\z-a ii \<R' i (A)}, (21) 
i=i 
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where C is the complex field and R'^A) = J2 I 

the deleted absolute row sums 

of A 

Since -R^(L) = ^ holds, (21) can be represented as follows: 

n 

G(L) = \j{zeC\ \z-£u\ <£ u }. (22) 
i=i 

Hence, we have 

Proposition 9. (1) The real part of every eigenvalue of L is nonnegative: every Gersgorin disc 
belongs to the right coordinate half-plane; 

(2) the intersection of all Gersgorin discs contains zero; 

(3) G(L) = \(z+ 1) max ta \z\<l\. 

y l<i<n > 

Obviously, the intersection of all Gersgorin discs consists of zero iff the digraph contain an 
undominated vertex. 

Consider the characteristic polynomial of L: 

i=0 

where Ei(L) is the sum of all principal minors of order i. By Theorem 6 in [2], Ei(L) = a ri 
for every i = 1, . . . ,n. Since every principal minor of order greater than n — v is zero, we have 



n—v . n—v 



p L (\) = (-i) n ~ v x v E a* {-X) n ~ v ~ l = x v E (-i)v, a— -\ 

i=0 i=0 
n—v 

Proposition 10. p' L {\) = A J2 °~ n -v-i ( — ^ s an annihilating polynomial for L. 

i=0 



8. ACCESSIBILITY VIA FORESTS AND DENSE FORESTS IN DIGRAPHS 

8.1. Forest accessibility 

The entries of Q(t) measure the proximity of the vertices of an undirected multigraph [14,15]. 
The matrix J = \im T ^ 00 Q T (T) was analyzed in [13] as the matrix of limiting accessibilities of a 
multidigraph. Here, we study the matrix Pi(t) = Q t {t) with r > as an accessibility measure 
for digraph vertices. By Theorem 4, the (i, j)-entry of this matrix is the total weight of out forests 
that "connect" i with j in the digraph where the weights of all arcs are multiplied by r. Along with 
P 1 (t), we consider the matrix of in forests P 2 (t). Its (i,j')-entry is the total weight of in forests 
(of the modified digraph) where j is a sink and i belongs to a tree converging to j. 

The following definition is formulated for an arbitrary vertex accessibility measure (formally, 
every square matrix of order n or, more precisely, the corresponding matrix-valued function of a 
digraph can be considered as such a measure). A measure P 2 is called to be dual to a measure P l 
if under the reversal of all arcs in an arbitrary digraph (provided that the weights of the arcs are 
preserved), the matrix of P 2 for the modified digraph coincides with P^[ calculated for the initial 
digraph. It follows from this definition that P 2 is dual to P 1 if and only if P 1 is dual to P 2 . In [15], 
three self-dual accessibility measures were studied. 

Let us check the satisfaction of the characteristic conditions listed below for P\{t) and P 2 {t). 
Triangle inequality for accessibility measures requires the symmetry of the corresponding matrix 
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(see, e.g., [26]). For that reason, we will check this condition for P 3 (r) = {P±(t) + P2(t) + Pii T ) + 
P 2 T (r))/4. 

Nonnegativity. For any digraph P, p^ > 0, i,j G V(r). 

Diagonal maximality. For any digraph P and any distinct i,j G V(P), 

(1) Pa > Pij and 

(2) Pll > Pji hold. 

Disconnection condition. For any digraph P and any i, j G V(P), p^- = if and only if j 
is unreachable from i. 

Triangle inequality for accessibility measures. For any digraph P and any i,j, k G V(r), 
Pij ^Pik ~Pjk — Pa holds. If, in addition, j = k and i ^ j, then the inequality is strict. 

Transit property. For any digraph P and any i,k,t G V(P), if P includes a path from i to 
k, i k t, and every path from i to t contains k, then (1) p ik > p it ; (2) p kt > p it . 

Monotonicity. Suppose that the weight of some arc e p kt in a digraph r increases. Then: 

(1) Ap kt > and for any i,j G V(r), / (k,t) implies Ap kt > Ap^-; 

(2) For any z G ^(i -1 ), if there is a path from k to £, and each path from k to i includes i, then 
(a) Ap fet > A Pkl and (6) A ftl > Ap ti ; 

(3) For any % G V(r'), if there is a path from i to k and every path from i to i includes k, then 
(a) Ap fet > A Pit and (6) Ap lt > Ap ik . 

The results of testing P 1 (t), P 2 (t), and P 3 (r) are collected in the following proposition. 

Proposition 11. The measures P\(t) and P2C7") are dua/ to each other for every r > 0. 
They satisfy nonnegativity, reversal property, disconnection condition, the first part of item 1, and 
item 2 of monotonicity. Moreover, P^r) satisfies items 1 of diagonal maximality and transit 
property; P 2 (t) satisfies items 2 0/ these conditions. With respect to the remaining statements of 
monotonicity, P\{t) satisfies items 2 and 3b, whereas P2(t) satisfies items 3 and 2b, and they both 
violate the second part of item 1. Furthermore, P\{t) breaks item 3a of monotonicity and item 2 
of transit property, whereas P 2 (r) breaks item 2a of monotonicity and item 1 of transit property. 
Triangle inequality for P 3 (r) is not satisfied. 

As was noted in [13], the limiting accessibility P = J T of a digraph does not completely corre- 
spond to the general concept of proximity. Notice that disconnection condition, which is satisfied 
for the limiting accessibility in one side only, is completely fulfilled for P^r) and P 2 (r). Moreover, 
Pi(t) and P 2 (r) obey a number of conditions which are satisfied by the limiting accessibility in 
the nonstrict form only. 3 

8.2. Accessibility via dense forests 

Now we consider a measure which is intermediate between the limiting accessibility (which 
depends on Q n -v only) and the forest accessibility Q{t) (which is a weighted sum of all matrices Qk). 
This new measure is determined by the matrices Q n - v -\ and Q n -v (or, equivalently, by the matrices 

3 By the nonstrict form of a condition we mean the result of substituting nonstrict inequalities (> and <) 
for the strict ones (> and <) in it. 
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J n -v-\ and Jn-v = J, which also determine L* as stated in Proposition 8). This measure can be 
also obtained by the inversion of L + a J with some values of a. 

Thus, consider the matrices R(a) = (r^) = (L + a J)' 1 with a > 0. Using Theorem 11 and 
Proposition 8, we have 

(L + a J)" 1 = L* + a- 1 J = ^zlj^^ + (V 1 - ^zl\ J. (23 ) 
If < a < an - v then, by (23), (L + a J) -1 is the sum of Q n - v -i and Q n -u with positive 

— v — 1 

coefficients. Spanning rooted forests (of an undirected multigraph) with n — v or n — v — 1 arcs are 
called in [15] dense forests, and the undirected counterpart of the accessibility measure (23) with 
< a < J 7 "-" - is called accessibility via dense forests. 

Consider two accessibility measures for digraphs: Pi(a) = R T (a), accessibility via dense diverg- 
ing forests and P 2 (a), accessibility via dense converging forests. 

An important property of the set of dense diverging forests is as follows. 

Proposition 12. 1. For any vertex i G V(r), there exists an out forest in T n -v-\ where i is 
a root. 2. For any path (chain subgraph) in r, there exists an out forest in T n -v-\ U Tn-v that 
contains this path. 

A similar proposition is true for converging forests. At the same time, the set of maximum out 
forests Tn-v and the set of maximum in forests do not have this property. For example, on Fig. 1 
in [13], no maximum out forest contains arc (4,2). 

We now test -Pi (a) and P 2 {ot). Similar to the previous consideration, triangle inequality for 
accessibility measures will be checked for the index P 3 (a) = (P 1 (a) + P[(a) + -P2( a ) + ^2 ( a ))/4> 
since this inequality requires the symmetry of the corresponding matrix. 

Proposition 13. For any a g] 0, a n - v / a n - v -\[, the measures P\(oi) and P2(&) are dual to each 
other. They satisfy nonnegativity and disconnection condition. Moreover, the nonstrict versions of 
items 1 of diagonal maximality and transit property are satisfied by i-* 1 (a), and items 2 of these 
conditions by P 2 {a). Both measures violate monotonicity. Triangle inequality for accessibility 
measures is not true for P 3 (a). 

CONCLUSION 

The normalized matrices of out forests are stochastic and determine the transition probabilities 
in the geometric observation model applied to the Markov chains related to the digraph under 
consideration. Various matrices of forests can be represented by simple polynomials of the Kirchhoff 
matrix. The Moore-Penrose generalized inverse L + and the group inverse L* of the Kirchhoff 
matrix L can be explicitly represented via L and the normalized matrix J of digraph's maximum 
out forests. The matrices of diverging and converging forests characterize the pairwise accessibility 
of vertices. These and other results enable one to consider the matrices of spanning forests as a 
useful tool for the analysis of digraph's structure. 

APPENDIX 

Proof of Proposition 1. 1. At first, let k < k' < n. r is constructed as follows. We draw 
a diverging star rooted at the first vertex and having k! — k leaf vertices. Also, we draw a path 
diverging from the root of the star and containing, in addition to the root, n — k' > 1 vertices that 
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are not included in the star. The remaining k — 1 vertices are left isolated. Then v = l + (k — 1) = k 
and v' = 1 + (A/ — k) + {k — 1) = k' , as required. 

If k' < k < n, then we draw a star converging to the first vertex and having k — k! other vertices. 
Also, we draw a path diverging to the center of the star and containing, in addition to the sink, 
n — k > 1 vertices that are not included in the star. The remaining k! — 1 vertices are left isolated. 
As well as in the first case, we have v = k and v' = k! . The second statement is obvious. 

Proof of Theorem 6. Since the spectral radius of P is 1, we have 

00 k — 1 

E((i-«)p) =(i-(i- q )py . 

k=0 

Using the formula of total probability and equations (3)-(5), we obtain 

P(a, g) = J2 P(k)P k = E ?(1 - Q) k P k = ?(/-(!- q)P) ~ (24) 

k=0 k=0 

= q(l -{l-q){I- aL)) ^ = q(ql + (1 - q)aL) ' 
= {l + ^^ L y = (I + tL)-i = Q(t). 

Proof of Proposition 3. In [13], we used the notion of weighted 2-digraph: it is a multidigraph 
with arc multiplicities no more than two. The weight of a 2-digraph is the product of the weights 
of its arcs. For a weighted digraph H and its vertices u,w G V(H), by H + (u,w) we denote 
the 2-digraph with vertex set V(H) and the arc multiset obtained from E(H) by the increment 
of the multiplicity of (u,w) by 1. Similarly, if H is a 2-digraph and u,w <G V(H), then by H' = 
H—(u, w) we denote the 2-digraph that differs from H in the multiplicity of (u, w) only: n'((u, w)) = 
max(n((«,w)) — 1,0). 

Now introduce the following notation. Let fl~^ s +(£, i) = {Fl~^ s + | Pl~^ S G !Fl~^ s }- For 
all i / j, by T 3 ^ s denote the set of out forests with k arcs where i is a root and s belongs to a tree 

diverging from j. Obviously, in such digraphs, s is unreachable from i whenever i ^ j. By 

we denote the set of all forests with k arcs where i is not a root, and s belongs to a tree diverging 

from j. These definitions induce the matrices {q^j) and (c/^j) with elements 

a k ., = e(T j ^ s ) 

<fij,i = e(H7')> s,j = l,...,n. 

Denote by cfi. the total weight of all diverging forests where i is unreachable from j. For all 
i,j€V(r), gf/=g& + 4- holds. 

Lemma 3. For any weighted digraph and all i, j = 1, . . . , n, we have: 

(1) Ifsjtp,i^j, and ( 8 ,i),(p,i)€E(r), then (^ k J s +(s, i)) f| {T^ p +(p, i)) = 0; 

(2) U (^CT+(M)) = ^+i- 
( s ,i)eE(r) 

Proof of Lemma 3. The first statement follows from the definition of out forest. Let us prove 
the second statement. Suppose that (s,i) G E(r) and JF^ S / 0. Then s is unreachable from i in 
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every forest Fk G J^IT* '■ After the addition of (s,i) to Fk, we obtain a diverging forest with k + 1 

arcs where i is reachable from j, i.e., Ffc + G 

Suppose now that F^+i G Fl+i- Let (P'*) ^ e ^ e unique arc directed to i in F^+i. Since 

i^fe+i - (P, *) G ^' ) 7 P , we have F fc+1 G i) C U + *))■ 

(«,0es(r) 

Corollary from Lemma 3. For any digraph and alH, j = 1, . . . , n and fc = 0, . . . ,n — u, 

= = < U (^',7* +(».<)))= E 4^7 s +(m)) = E 

(s,i)eB(r) (s,j)es(r) ( s ,j)e£(r) 

We now continue proving Proposition 3. Let LQk = (a^ ). Then for alH = 1, . . . , n, 

n n 

4=^2^ = J2 e ^si+^Qu = E^ - E^^ = S e «(«« - ( 25 ) 

s=l s^i s^i s^i s=l 

n n 

= E e «(«a + ^ " 9*) = E e «4 = ^fc+i > 0- 

The statement of Proposition 3 with respect to the diagonal entries of LQk is proved. 
For j 7^ i we have 

n n n 

a jj = ^ ' ^isQsj = y ' ^isQsj ~ y ' ^isQij ^ ] ^siQij ~ ^ ] ^siQsj 
s=l s^i Sy^i s=l s=l 



= ]T £si g§ - ]T £«(<&i + = E e »«?j " E £ «<i " E £ si1sj,i 

s=l s=l s=l s=l s=l 

n 

= ^i)-e(^)-E e -4" ( 26 ) 

s=l 

where e(Qi) = Y^=i e siQij, z^Q-i) = Y11=\ £ siq k s j- v Gi is the multiset of 2-subgraphs obtained by the 
addition of all arcs (s, i) G E(r) to all possible forests in and Q2 is the multiset of 2-subgraphs 

obtained by the addition of all arcs (s, i) G E(r) of r to all possible forests in J-V-^ i '■ The multiset 
£1 consists of the pairs (H,m(H)) where m(H) > 1 is the multiplicity of H in Q\\ Q2 consists of 
the pairs (H,ri2(H)) where ri2(H) is the multiplicity of H in Q 2 - The notation H G Q\ means here 
ni (H) > 0. 

Prove that Q\ = Q2- in every H G Qi, two arcs are directed to i: (s,i) and (p, i) such that 
(p, i) G E(F^ 1 ), since j 7^ i. Consider the digraph H — (p,i). It is a forest where i is not a root 
and j is the root of the tree that contains p. Then H — (p, i) G J- J k ^ P and H G T 3 k ^ v +(Pj *) ^ £?2- 
Conversely, suppose that = F^ s + (s,i) G C/2- Consider — (p,i) such that (p, i) G E(FjT^ s ). 

We have — (p,i) G J 7 ^ 1 , hence, H = F^ 1 + (p, i) G C/i. For the multisets Q 1 and C?2 , the 
multiplicities of all elements, n\{H) and n,2(H), do not exceeded two. Otherwise, H — (s,i) would 
not be a forest. Prove that for any H, n\(H) = 2 iff ri2(H) = 2. Suppose that H G Q\ and 
ni(H) = 2. Then there are vertices si,S2 G V(r) such that si 7^ S2 5 E(H) contains and 
(s2,i), and {i? — (si,i), H — (s2,i)} Q T 3 ^ 1 ■ Since s\ is reachable from j in H — (s2,i), si is 
reachable from j in H too. Therefore, si is reachable from j in H — (s\,i) as well. Similarly, 
S2 is reachable from j in H — (s2,i)- We obtain H — (s\,i) C T 3 ^ 1 and H — (s2,i) C T 3 k ~j S2 , 
consequently, ri2(H) = 2. 

Suppose that H G Q2 and ri2(H) = 2. Then for some distinct vertices si and S2, there exist 
G F^ S1 and FfV 52 G jr^ S2 such that F = Ff'r Sl + = FfV 52 + (s 2 ,i). Note that 
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(A) F^ S1 contains the arc (s2,i), therefore, S2 is not reachable from % in this forest. Then (B) S2 
is reachable from j in F^ S1 , since otherwise S2 would be reachable from j in H = F^ S1 + (si,i), 
which is false. By (A) and (B), F^U Sl e T{^\ and similarly, F^ 82 G JFjP*. Hence, m(H) = 2. 
By (26) and Corollary from Lemma 3, at j ^ i we have 

4 = -]>>^ = -^ +1 <0. (27) 

Proposition 3 is proved. 

Proof of Theorem 7. Let us represent (14) as follows: 

' Qo = cr Q I = I, 
Qi = (T 1 I — LQ , 



(28) 



Qk = <* k l - LQk-i, 

„ Qn—v = ® n—v I LQ n - v -\. 

By substituting each equation in the subsequent one, for all k = 2, . . . , n — v we have 

Q h = a h I-L(a k _ 1 I-...-L{a 1 I-L<T Q /)...) 

fc 

= a k I - a k _, L + ... + (-l) k a L k = £ (-L)\ 

i=0 

Theorem 7 is proved. 

Proof of Lemma 1. The ith row sum of LQ k = (a k j) is 

n n n n n n 

E 4 = E E = E ^ E A = E = °- 

j=l j=l s=l s=l j=l s=l 

Proof of Corollary 1 from Theorem 7. Multiplying both sides of (11) by any matrix 
that commutes with L and using distributivity and associativity of matrix operations, we get the 
required statement. 

Proof of Corollary 2 from Theorem 7. Consider (25) at k = n — v. By virtue of Proposi- 
tion 2, for every (s,i) G E(r), e s iq^~ v = holds, i.e., a™~ v = for all i € {1, . . . , n}. In this way, 
Corollary 2 is derived from a^~ v = 0, inequality (27), Lemma 1, and Corollary 1 from Theorem 7. 

Proof of Corollary 3 from Theorem 7. Postmultiplying both sides of (11) by J and using 
LJ = (Corollary 2 from Theorem 7) provides Q k J = a k J. By commutativity, JJ k = JkJ = J 
holds. Using Theorem 4, we also get Q(t)J = JQ(t) = J for any r > 0. 

Proof of Theorem 8. Substituting (11) in (15) gives 

n—v k n—v n—v 

s (j + l)- 1 = y: E ^(-ly = E E °k- t (-Ly 

k=0 i=0 i=0 k=i 

n—v n—v—i n—v 

= E E "H^EwiW 

i=0 j=0 i=0 
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Proof of Corollary from Theorem 8. Observe that the digraph resulting from r by multi- 
plying the weights of all arcs by r has the Kirchhoff matrix tL, and its total weight of out forests 
with j arcs is t 3 o j . Hence, the required statement follows from Theorem 8. 

Proof of Theorem 10. We first prove the following lemma. 



Lemma 4. If a maximum out forest of a digraph r is a tree, 
r is 1, then L + J is nonsingular. 



i.e.. 



the out forest dimension of 



Proof of Lemma 4. Assume, on the contrary, that det(L + J) = 0. Then there exists a vector 
b = (&!,... , b n ) T + such that (L + J) T b = 0, where = (0, . . . , 0) T . 

Since, for every digraph of out forest dimension 1, every column of J consists of equal entries 
(item 3 of Theorem 1), we obtain 



(L + Jfb 



Mil + • • • + b n i n i + bi Jn + . . . + b n Jni 
M12 + • • • + b n i n2 + hJ 12 + . . . + b n J n2 

b\l\ n + • • • + b n l nn + b\J\ n + . . . + b n J nn 



Mn + • 


• + b n £ nl + Jn(6i + . 


■ + b n ) 







Ml2 + • 


• + Wn2 + Jl2{h + . 


■ + b n ) 







Mln + • 


• + b n i nn + Jinih + . 


■ + b n ) 








(29) 



Adding up the components of the vectors that form the last equality and using the identities 

n _ n 

Jik = 1 (item 1 of Theorem 1) and J2^ik = 0, i = 1, . . . , n, we deduce b\ + . . . + b n = 0. 

k=l k=l 

Replacing the last equation of (29) with this equality, we obtain the following system of equations 
in &i,...,&„: 

Mn + • • • + b n e nl = o 

M12 + • • • + b n £ n2 = 

(30) 

Mln-l + • • • + bnlnn-! = 
Ui + . . . + b n =0. 

Let M be the transposed matrix of coefficients of the system (30). Expand the determinant 
of M in the last column (which consists of ones). By the matrix-tree theorem for digraphs (see, 
e.g., [1, 16]), if there exists a spanning tree diverging from i, then the cofactor of £{ n in L — which 
equals the cofactor of the ith entry of the last column of M — is positive, whereas, in the opposite 
case, it is zero. By the hypothesis of this lemma, the maximum out forests of r are diverging 
trees. Therefore, the cofactor of at least one entry of the last column of M is positive. Hence, 
det M expanded as above is also positive. Thereby, rankM = n, whence the unique solution of the 
system (30) is b\ = ■ ■ ■ = b n = 0. This contradiction proves Lemma 4. 

Suppose that V\.,...,V P are the vertex sets of the weak components of r. Without loss of 
generality, we assume that the vertices of V\ are numbered first, the vertices of V 2 are numbered 
next, etc. (at any other numeration, the corresponding permutation of the rows and columns 
preserves the rank of L + J). By virtue of Theorem 2' in [13], L + J is a block diagonal matrix 
with p blocks, which is expressed in the following manner [25]: 



L + J 



(31) 
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where submatrix A ss corresponds to the sth weak component. A block diagonal matrix L + J is 
nonsingular if and only if every A ss is nonsingular, moreover, 

p 

det(L + J) = Yl detA ss . (32) 

s=l 

Let r s be the restriction of r to V s , s G {1, . . . L^, Q^ s 2„, and J 1 ^ will be the matrices 
constructed for r s . Then 

A ss = + J (s) . (33) 

Indeed, l/ s ) coincides with the sth block of L, whereas, by item 5 of Theorem 1, the sth block 
of Q n - v is proportional to Qnl v (the proportionality factor being the total weight of all forests in 
r^Vs)- Consequently, the sth block of J coincides with J^ s \ 

By virtue of (32) and (33), it is sufficient to prove the statement of the theorem for the case of 
p = 1. Let r consist of a single weak component. Suppose, without loss of generality, that V(r) 
is indexed in such a way that the first numbers are attached to the vertices in K\ , the subsequent 
numbers to the vertices in K 2 , etc. The last numbers are given to the vertices in V(r)\K. By 
item 2 of Theorem 1, L + J is a block lower triangular matrix with v + 1 blocks: K\ corresponds 
to the first block, K v corresponds to the vth. block; V(r)\K corresponds to the last block. 

The determinant of the block triangular matrix L + J is the product of the determinants of 
its diagonal blocks, whereas its rank is no less than the sum of ranks of the diagonal blocks (see, 
e.g., [25]). Note that, as well as in the case of weak components of r (see above), the block of L + J 
corresponding to an undominated knot Ki, coincides with the matrix L(i~^ i ) + J(i^ i ) constructed 
for r^i- To demonstrate this, it suffices to use item 3 of Theorem 1. Every maximum out forest of 
an undominated knot is a diverging tree (out arborescence). Therefore, the nonsingularity of the 
diagonal blocks of L + J that correspond to the undominated knots follows from Lemma 4. 

The (v + l)st block of L + J coincides with the corresponding block of L, since the last block of 
J is zero (by item 2 of Theorem 1). The last block of L is nonsingular. Indeed, by Theorem 6 in [2], 
its determinant is the weight of the set of out forests in r where K is the set of roots. This weight 
is strictly positive, since the indicated set of forests is nonempty: such forests can be obtained from 
the maximum out forests of r by the removal of all arcs between the vertices within K. 

Thus, the diagonal blocks of L + J are nonsingular, hence, L + J is nonsingular. The theorem 
is proved. 

Proof of Corollary from Theorem 10. First, we prove the following lemma. 

Lemma 5. For every weighted digraph r of out forest dimension 1 and any a ^ 0, the matrix 
L + a J is nonsingular. 

Proof of Lemma 5. This lemma is proved by the same argument as Lemma 4 with the only 
difference that the analogue of (29) takes here a more general form: 





h£n + . 


■ + Mni + atJu(bi + . 


■ + b n ) 







(L + a J) T h = 


W12 + • 


• + b n £ n2 + aJ 12 (b! + . 


■ + b n ) 









hhn + ■ 


• + Wnn + a Ji n (61 + . 


■ ■ + b n ) 








To complete the proof of Corollary from Theorem 10, note that for any a 7^ 0, the matrix 
L + aJ, as well as L + J, is a block lower triangular matrix with v + 1 blocks. By item 2 of 
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Theorem 1, its (v + l)st diagonal block coincides with the corresponding block of L + J. Using 
Lemma 5, we conclude that the other diagonal blocks are also nonsingular. 

Proof of Theorem 11. By Theorems 2 and 3, (L + J) J = J. Premultiplying both sides of 
this identity by (L + Jy 1 (which exists by Theorem 10), we have 

J=(L + jy 1 J. (34) 

Similarly, 

J{L + J)- 1 = J . (35) 
Denote Q = (L + J)" 1 — J . Using (34) and Theorems 2 and 3, we obtain 

QL = (L + J)~ 1 L — J L = (L + J)~ 1 {L + J - J) = I - (L + J)' 1 J = I - J . (36) 

Similarly, 

LQ = I-J. (37) 

By (34) and Theorem 2, 

QJ= (L + J)' 1 J- J 2 = 0. 

Consequently, for any a^O, 

(Q + a' 1 J)(L + a J) = I - J + J = I. 
Hence, Q + a -1 J = (L + a J) -1 and 

Q = (L + a J)' 1 - a' 1 J (38) 

for any 

By (36) and (37), LQ = QL = I — J, thus Q = (L + a J) -1 — a" 1 J satisfies condition (5) in the 
definition of group inverse. Let us prove that conditions (1) and (2) (common with the definition of 
Moore-Penrose inverse) are also fulfilled. Making use of Theorems 2 and 3 and identities (35)-(38), 
we obtain 

LQL = L(I - J) = L, 

QLQ = (I - J)Q = Q - JQ = Q - J (L + Jy 1 + J 2 = Q - J + J = Q, 

which completes the proof. 

Proof of Proposition 7. Using Theorems 2, 3, and 5 and identities (I + rLy 1 ] = J 
(Corollary 3 from Theorem 7) and (34), we obtain 

( Km r ((/ + riy 1 - J) + J) (L + J) 

= lim r ((/ + rLy 1 L + (/ + riy 1 J-JL-J 2 ) + JL + J 2 
= lim t(I + rLy 1 L + J = lim (/ + rL) _1 (7 + tL-I) + J 

T— »0O T— »00 

= /- lim (I + rLy 1 + J = 1. 

Postmultiplying the first and last expressions by (L + Jy 1 and using Theorem 11 we obtain the 
required equation. 
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Proof of Proposition 8. Using Proposition 7 and Theorem 4, for any i,j G V(r), find the 
limit 



£ ti = ^r^^r) - ) = Um r 



fc=0 



E r k a k 



Jij 



n—v 

. E 

\ k=0 

n—v . . n—v 

E T^e(FD - E r k ^a k J y 

v fe=0 fc=0 

nm 

E r*<7 fc 
fc=0 

n— 1 . . . n—v—1 

E ^(O + ^MO- E T^^Jy-T-^V^Jy 
fc=0 fc=0 



= lim 



E T k (J k 
k=0 



By the definition of J, we have r™ — t" v+1 a n _ v Jij = 0. Therefore, in view of 

(J n - V / 0, this yields 

pf _ £ (-7"" n-v-l) ~ ff n-«-l ^jj 
n—v 

which completes the proof of Proposition 8. 

Proof of Lemma 2. By virtue of the identity J L = (Theorem 3), matrices LL T and ZZ T = 
J T J+LL T commute, i.e., LL T (J T J+LL T ) = (J T J+LL T )LL T = (LL T ) 2 . Premultiplying and 
postmultiplying both sides of the first equality by (ZZ T )~ 1 = (J T J +LL T )~ 1 , we obtain the desired 
(J J +LL T )~ 1 LL T = LL T (J J +LL T )~ 1 . The second statement is proved similarly. 

Proof of Theorem 12. Let = L T {Z Z T )~ l . Prove that the following four conditions are 
satisfied: 

(1) LL^L = L; 

(2) L< + >LL<+> = L<+>; 

(3) LL (+) is symmetric; 

(4) L (+) L is symmetric. 

Condition 1. Using Lemma 2 and the identity J L = (Theorem 3), we obtain 

LL M L = LL T {J T J+LL T )- l L = {{LL T + J T J) - J T J)(J T J +LL T )~ 1 L 
= (I- J T J{J r J+LL T )- 1 )L = L — {J T J +LL T )- 1 J T JL = L. 

Condition 2. Using Lemma 2 and the identity L T J T = 0, we have 

L (+) LL <+) = L {+) LL T (J T J +LL T )~ 1 = L<+>((LL T + J T J) - J T J){J T J +LL T y 1 

= L (+) (I - J T J{J r J+LL T )- 1 ) = L (+) - L T (J T J+LL T )- 1 J T J(J T J +LL T y l 
= L w - L T J T J{J T J+LL T )- 2 = L (+ \ 

Condition 3. By Corollary from Lemma 2, the matrix LL (+) = LL T (J T J +LL T )~ 1 is symmetric. 
Condition 4. Since (J T J +LL T )~ 1 is symmetric, L (+) L = L T (J T J +LL T )~ 1 L is symmetric as 
well (see, e.g., [25], Theorem 4.1.3). 
Theorem 12 is proved. 
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n—v 

Proof of Proposition 10. By Theorem 7, Q n -v = °~ n -v-i (~ • Using Corollary 2 from 

i=0 

Theorem 7, we obtain 

n—v 



p' L (L) = LJ2 °n-v-i (-LY = LQ n - v = 0. 
i=0 

Proof of Proposition 11. Under the reversal of all arcs in a digraph, diverging and converging 
forests change places. Therefore, Pi(t) and -P 2 ( r ) are dual to each other. 

Nonnegativity and disconnection condition for Pi(t) and -P2( r ) and also item 1 of diagonal 
maximality for -P 1 (r) and item 2 of diagonal maximality for P 2 ij) are proved with the help of 
Theorem 4 by the same argument as the corresponding conditions in [14]. 

For any i,k, t € V(G), if r contains a path from % to k, i ^ k ^ t, and every path from i to 
t includes k, then J^^ 1 C j^^ k . This inclusion and Theorem 4 imply the fulfillment of item 1 of 
transit property for -P 1 (r). Item 2 of transit property for P 2 ( r ) ^ s P r oved similarly. 

It was established in [13] (in the proof of Proposition 17) that if the weight of some arc (k, t) is 
increased by Ae^t and the weights of all other arcs are preserved, then the increments of the entries 
of the matrix Pi(t) = (Pif ( T )j are expressed as follows: 

Ap/(r) = — — ^y— ^y— , i,j = l,...,n. (39) 

By reversal property, the corresponding formula for -P 2 ( r ) = {pT^ ( T )) ls 

A Pij( T ) = ; m t,j = i,...,n, (40) 

where pf^ {t) is the total weight of in forests wherein j belongs to a tree converging to i. 



Transcribe (39) for Ap^(r): 



A^(r) 



«^ J 1 ) 



Since PjP(r) > (r) for all i, j = 1, . . . , n, we deduce Ap^ (t) > 0, i.e., P 1 (t) satisfies the first 
part of item 1 of monotonicity. 

Now represent A p^ (r) using (40) : 

( 2 ), , ptl{T)(p${T)-p${T) 



^Pl7(r) 



Since p^(r) > p^(r) for all i, j = 1, . . . ,n, we have Ap[ 2 ( ' (r) > 0, i.e., P 2 (t) also satisfies the 
first part of item 1 of monotonicity. 

Transcribing (39) for Ap^ (r), Ap^(r), and Ap|-'(r) and using item 1 of diagonal maximality, 
we conclude that P-y{r) satisfies item 2 of monotonicity. Comparing the expressions for Ap^(r) 
and Ap-P(r) and using item 1 of transit property, we obtain item 3b of monotonicity. Along the 
same lines, the statement of item 3a holds if and only if (p^ (t) —p^t ( T )) ~~ {pik ( T ) ~~ Pit ( T )) > 
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Since triangle inequality for accessibility measures is violated by P 1 (r) (as well as by all asymmetric 
measures), item 3a of monotonicity is not true. Similarly, P 2 (r) satisfies items 3 and 26, but violates 
item 2a of monotonicity. 

Consider the digraph with vertex set {i,j,k,t}, arc set {(i,k),(k,t),(t,j)}, and arc weights 
e(i, k) = 4, e(k,t) = 1, and e(t,j) = 4. Then P x (l) is given by 



^i(l) 



i j k t 

1 0.32 0.8 0.4 
0.2 
0.08 0.2 0.1 
0.4 0.5 



(41) 



Since p\l\l) > p ( kt '(l), item 2 of transit property is not satisfied. 
Compare the increments Apj^(l) and Ap^\l) for an arbitrary Ae kt > 0: 



,(!) 



P ( t ]\l)(p^(l)-P ( u ) (l] 



Since pg)(l) -pg^l)) = 0.5(0.2 



0.1) 



0.05, pg^l) (p. 



(i) ( 



0.4(0.? 



0. 4) = 0.16 and the common denominator is positive, Apj^(l) < Ap^\l) follows, i.e., the second 
part of item 1 of monotonicity is not satisfied. 

It is easy to verify that item 1 of transit property and the second part of item 1 of monotonicity 
are violated by P 2 (r) m the same example. 

Let us show that P^(t) does not satisfy triangle inequality for accessibility measures. Consider 
the digraph with vertex set {i,j,k,t}, arc set {(i,j),(J,k), (k,t),(t,i)}, and arc weights s(i,j) = 

1, e(j,k) = 10, e(k,t) = 10, and e(t,i) = 1. Here, P 3 (l) is as follows: 



k 



0.6302 0.2233 0.1693 0.2233 
0.2233 0.3724 0.1823 0.2747 
0.1693 0.1823 0.1146 0.1823 
0.2233 0.2747 0.1823 0.3724 



(42) 



Triangle inequality for accessibility measures is violated, because p ki (I) + p k j (1) — p\j (1) > 

(3) 

p kk (1). Symmetry is a necessary condition of triangle inequality for accessibility measures. There- 
fore, Pi(t) and P 2 (t) do not satisfy this inequality either. 

Proof of Proposition 12. 1. To construct the required forest, it suffices to take a maximum 
out forest and to delete any arc directed to % (if such an arc exists) or any arc in the opposite case. 

2. For the given path (chain subgraph) and any maximum out forest in P, consider their join 
and remove all arcs of the out forest that are directed to the vertices of the path but do not belong 
to the path. The resulting subgraph contains neither circuits nor vertices with indegree greater 
than one, i.e., it is an out forest. It contains at least n — v — 1 arcs, consequently it belongs to 

J~n—v—l U J-n—v ■ 
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Proof of Proposition 13. Under the reversal of all arcs in P, diverging and converging forests 
change places. Hence, Pi (a) and P 2 (a) are dual to each other. Disconnection condition follows 
from item 2 of Proposition 12. 

Nonnegativity follows from the fact that the entries of J n - V and J n - v -\ are proportional to the 
weights of some sets and from the nonnegativity of the coefficients in (23). 

Item 1 of diagonal maximality and item 1 of transit property for Pi(a) in the nonstrict form 
(as well as items 2 of these conditions for P 2 (a) i n the nonstrict form) follow from the nonstrict 
inclusion of the sets of forests that determine the entries of matrices Q n - V and Q n - v -i under 
comparison. 

Diagonal maximality (in the strict form) is violated, for example, for digraphs with vertex set 
{j,i,k,t}, arc set {(j,i), (i, k)(k,t)}, and arc weights e(j,i) = 4, e(i,k) = 1, and e(k,t) = 1. The 
matrices of out forests Q n -v an d Q n -v-i, the matrices of in forests S n - V and S n - v -i, and the 
matrices Pi (a) and P 2 (a) at a = 4/13 are as follows: 



Qn- 



j i kt 




j i kt 




400 


3 


900 


3 


400 


i Qn-v-l = 


8 10 


i 


400 


k 


4 14 


k 


400 


t 


14 4 


t 



Sn—v — 



j i kt 




00 4 


3 


00 4 


i 


00 4 


k 


00 4 


t 



Sn—v—l 



j i kt 

1440 
044 1 
04 5 
00 9 



3 


i k t 




3.25 3 2 1 







0.25 0.25 0.25 





1 1 







1 






j i kt 






0.25 1 1 1 


3 




1 1 1.25 


i 




1 2.25 


k 




3.25 


t 



In this example, item 1 of transit property for Pi(a) is violated (since p^ (a) = p^\a)), and 
so is item 2 of transit property for P 2 (a) (since p^(a) = p^(a)). 

Let us demonstrate that triangle inequality for accessibility measures is not satisfied by P^ia) = 
(Pi (a) + P[(a) + P 2 (a) + P 2 T (a))/4 in this example. Indeed, 



P 3 (a) 



j i k t 

1.75 1 0.75 0.5 
1 0.625 0.3125 0.375 
0.75 0.3125 1 0.8125 
0.5 0.375 0.8125 2.125 



and since p^{a) + p^\a) — pff (a) = 0.875 > p\f' '(a) = 0.625, this condition is violated. 

Monotonicity is not satisfied on undirected graphs (see Proposition 10 in [15]), hence, it is 
violated for digraphs also. 



REFERENCES 

1. Harary, F., Norman, R.Z., and Cartwright, D., Structural Models: An Introduction to the Theory of 
Directed Graphs, New York: Wiley, 1965. 



AUTOMATION AND REMOTE CONTROL Vol. 62 No. 3 2001 



466 AGAEV, CHEBOTAREV 

2. Fiedler, M. and Sedlacek, J., O IV-basich Orientovanych Grafu, Casopis Pest. Mat., 1958, vol. 83, 
pp. 214-225. 

3. Kelmans, A.K. and Chelnokov, V.M., A Certain Polynomial of a Graph and Graphs with an Extremal 
Number of Trees, J. Comb. Theory, 1974, scr. B, vol. 16, pp. 197-214. 

4. Liu, C.J. and Chow, Yu., Enumeration of Forests in a Graph, Proc. Amer. Math. Soc, 1981, vol. 83, 
pp. 659-663. 

5. Chaiken, S., A Combinatorial Proof of the All Minors Matrix Tree Theorem, SIAM J. Discr. Math., 
1982, vol. 3, pp. 319-329. 

6. Myrvold, W., Counting fc-component Forests of a Graph, Networks, 1992, vol. 22, pp. 647-652. 

7. Bapat, R.B. and Constantino, G., An Enumerating Function for Spanning Forests with Color Restric- 
tions, Linear Algebra AppL, 1992, vol. 173, pp. 231-237. 

8. Takacs, L., Enumeration of Rooted Trees and Forests, Math. Scientist., 1993, vol. 18, pp. 1-10. 

9. Erdos, P.L., A New Bijection on Rooted Forests, Discrete Math., 1993, vol. Ill, pp. 179-188. 

10. Merris, R., Doubly Stochastic Graph Matrices, Univ. Beograd. Publ. Elektrotehn. Fak. (Ser. Mat.), 1997, 
vol. 8, pp. 64-71. 

11. Merris, R., Doubly Stochastic Graph Matrices II, Lin. Multilin. Algebra, 1998, vol. 45, pp. 275-285. 

12. Chebotarev, P.Yu. and Agaev, R.P., The Matrix of Maximum Out Forests and Structural Properties 
of Systems Modeled by Digraphs, Modelling and Simulation of Systems, MOSIS-2000, 34th Spring Int. 
Con]., Ostrava, 2000, vol. 1, pp. 101-106. 

13. Agaev, R.P. and Chebotarev, P.Yu., The Matrix of Maximum Out Forests and Its Applications, Avtom. 
Telemekh., 2000, no. 9, pp. 15-43 [Automat. Remote Control, 2000, vol. 61, no. 9, pp. 1424-1450.] 

14. Chebotarev, P.Yu. and Shamis, E.V., The Matrix-Forest Theorem and Measuring Relations in Small 
Social Groups, Avtom. Telemekh., 1997, no. 9, pp. 124-136. [Automat. Remote Control, 1998, vol. 58, 
no. 9, pp. 1505-1514]. 

15. Chebotarev, P.Yu. and Shamis, E.V., On Proximity Measures for Graph Vertices, Avtom. Telemekh., 
1998, no. 10, pp. 113-133 [Automat. Remote Control, 1998, vol. 59, no. 10, pp. 1443-1459]. 

16. Harary, F., Graph Theory, Reading, Mass.: Addison- Wesley, 1969. 

17. Tutte, W.T., Graph Theory, Reading, Mass.: Addison- Wesley, 1984. 

18. Zykov, A. A., Teoriya konechnykh grafov (Theory of Finite Graphs), Novosibirsk: Nauka, 1969. 

19. Faddeev, D.K. and Faddceva, V.N., Computational Methods of Linear Algebra, San Francisco: 
W.H. Freeman and Co., 1963. 

20. Berman, A. and Plemmons, R., Nonnegative Matrices in the Mathematical Sciences, New York: Aca- 
demic Press, 1979. 

21. Gelfand, I.M., Lektsii po lineinoi algebre (Lectures on Linear Algebra), Moscow: Nauka, 1971. 

22. Campbell, S.L. and Meyer, CD., Generalized Inverses of Linear Transformations, London: Pitman, 
1979. 

23. Chebotarev, P.Yu. and Shamis, E.V., Preference Fusion when the Number of Alternatives Exceeds Two: 
Indirect Scoring Procedures, J. Franklin Inst., 1999, vol. 36, pp. 205-226. 

24. Hall, K., An r-dimensional Quadratic Placement Problem, Management Sci., 1970, vol. 17, pp. 219-229. 

25. Horn, R.A. and Johnson, C.R., Matrix Analysis, Cambridge: Cambridge Univ. Press, 1986. 

26. Chebotarev, P.Yu. and Shamis, E.V., On a Duality Between Metrics and S-proximities, Avtom. Tele- 
mekh., 1998, no. 4, pp. 204-209 [Automat. Remote Control, 1998, vol. 59, no. 4, pp. 608-612; Erratum, 
1998, vol. 59, no. 10, p. 1501]. 

This paper was recommended for publication by V.A. Lototskii, a member of the Editorial Board. 



AUTOMATION AND REMOTE CONTROL Vol. 62 No. 3 2001 



